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We consider models in which a dark-matter particle decays to a slightly less massive daughter 
particle and a noninteracting massless particle. The decay gives the daughter particle a small velocity 
kick. Self-gravitating dark-matter halos that have a virial velocity smaller than this velocity kick 
may be disrupted by these particle decays, while those with larger virial velocities will be heated. 
We use numerical simulations to follow the detailed evolution of the total mass and density profile 
of self-gravitating systems composed of particles that undergo such velocity kicks as a function of 
the kick speed (relative to the virial velocity) and the decay time (relative to the dynamical time). 
We show how these decays will affect the halo mass-concentration relation and mass function. 
Using measurements of the halo mass-concentration relation and galaxy-cluster mass function to 
constrain the lifetime-kick-velocity parameter space for decaying dark matter, we find roughly that 
the observations rule out the combination of kick velocities greater than 100 km s _1 and decay times 
less than a few times the age of the Universe. 
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I. INTRODUCTION 



There is a good consensus from many types of obser- 
vations that dark matter makes up ~ 25% of the mass- 
energy density of the Universe 1-5]. However, the na- 
ture of dark matter is unknown. The most popular class 
of candidates is the weakly interacting massive particle 
(WIMP), since such particles may be produced thermally 
in the early Universe at the right abundance and with be- 
havior consistent with a large set of cosmological observa- 
tions • The canonical WIMP is electrically neutral, 
and it is stable. Once produced in the early Universe, it 
interacts with itself and with ordinary matter only grav- 
itationally. If the WIMP is the dark matter, then the 
dark halos of galaxies and galaxy clusters are described 
by a collisionless gas of self-gravitating WIMPs. 

However, there is both theoretical and observational 
room for other types of dark-matter candidates. While 
observations are consistent with WIMPs, the observa- 
tions do not require the dark matter to be WIMPs. More- 
over, some observations — e.g., of the inner mass distribu- 
tion in galaxies and/or the abundance subhalos in dark- 
matter halos [ll[ — have inspired theoretical searches for 
dark-matter candidates with properties beyond those of 
the canonical WIMP. There is a large amount of lit- 
erature on dark-matter particles with additional physi- 
cal properties beyond those of the canonical stable colli- 
sionless WIMP. Examples include (but are by no means 



limited to) particles with small electric charges [12| or 
dipoles [lj|, self-interacting particles [Til, or particles 
with long-range forces |15l - fl9| (see Refs. 2CJ, [21| for re- 
cent reviews). 

In the spirit of these lines of investigation, we consider 
in this work a neutral dark-matter particle X of mass 
Mx that decays with lifetime r to a slightly less massive 
neutral particle Y of mass My = Mx(l — e) with e <C 
1 and an effectively massless particle £ which is itself 
assumed to be noninteracting. We imagine that such a 
scenario may arise in some implementations of models of 
inelastic dark matter 
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When the X particle decays, the daughter Y particle 
receives a nonrelativistic velocity kick ~ ec. Now sup- 
pose that these particles make up a self-gravitating halo 
of virial velocity i> v ; r . If <C v v ; r , then the halo may 
be heated slightly, and its mass distribution thus rear- 
ranged slightly, by the velocity kicks imparted to the Y 
particles once most of the X particles have decayed. If, 
on the other hand, 3> v V ir, then these halos will be 
completely disrupted after most of the X particles have 
decayed. There should thus be no halos with v V i r <C i>k 
if these particles decay with lifetimes r small compared 
with the age tjj of the Universe, and this has been pos- 
tulated as a possible explanation for the low number of 
dwarf galaxies relative to that expected from dissipation- 
less dark-matter simulations [23|, [24|. Alternatively, if 
fk Vvir but t 3> tjij then only a small fraction (~ t/<h) 
of the halo particles will be kicked out of the halo. The 
resulting halo mass and mass distribution may then be af- 
fected slightly without being completely disrupted. The 
canonical- WIMP halo is, of course, recovered in the lim- 
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its r — > oo and/or e — > 0. 

In this paper, we report on simulations of self- 
gravitating halos with decaying particles. While the re- 
sults of these decays can be understood in the limits 
r > £h and Vk v v i T [25[ using the adiabatic-expansion 
model [1^, |27j . detailed evolution of the halo over the 
full range of the r-i>k parameter space requires numeri- 
cal simulation. We use the simulations and observations 
of the galaxy-cluster mass function and the halo mass- 
concentration relation to constrain the decay parameter 
space. Our central results are presented in Fig. [31 which 
shows that (roughly) the combination of decay times less 
than a few times the age of the Universe and kick veloc- 
ities greater than 100 km s _1 are ruled out. 

The outline of the paper is as follows. In Sec. HH we 
describe our simulations. In Sec. IIII1 we characterize the 
evolution of dark-matter halos as a function of the de- 
cay parameters using the simulations, and show how the 
simulations, in conjunction with observational and theo- 
retical (in the context of WIMPs) determinations of the 
cluster mass function and the mass-concentration rela- 
tion, allow us to constrain the decay parameter space. In 
Sec. II VI we discuss various aspects of our findings, includ- 
ing resolution effects (with an eye towards the require- 
ments for cosmological simulations) and observational bi- 
ases. We summarize our work in Sec. [V] 



II. SIMULATIONS 

We simulate isolated equilibrium dark-matter halos us- 
ing the parallel TV-body code GADGET-2 in its New- 
tonian, noncosmological mode [HI, [29j]. We have mod- 
ified GADGET-2 in the following way to handle de- 
cays. The maximum time step Ai max allowed in the 
simulations is chosen such that the probability for a de- 
cay within that maximum time step is well-described by 
-fmax = Ai max /r < 1, At each time step, there is a 
Monte-Carlo simulation of decays, with the decay prob- 
ability being P = At/r for each particle. If a particle is 
designated for decay, it receives a kick speed «k in a ran- 
dom direction, and is flagged to ensure that the particle 
cannot decay again. Because we assume the mass differ- 
ence is negligible throughout the simulations (e < 10~ 3 ), 
we maintain the mass of the particles; the change in the 
kinetic energy of the particle due to the kick is much 
greater than the change to the potential energy due to 
the decrease in particle mass. The massless £ particle 
escapes from the halo and is not of interest in this calcu- 
lation. 

Our halos are spherically symmetric, and the density 
p(r) as a function of galactocentric radius r is taken to 
be the Navarro-Frenk- White (NFW) form (30l[3l|, 



p(r) = 



r r 
- 1 + - 



(1) 



such that the halo concentration is 



(2) 



where i? v ir is the virial radius, or size, of the halo. We 
define virial quantities with respect to the spherical top- 
hat overdensity, such that the virial mass 



4?r 



A vPm R 



3 

vir 3 



(3) 



where p m is the mean matter density of the Universe, and 
A v = (18tt 2 - 39a; - 82x 2 )/Q I[L (z) in a ACDM Universe, 
where x = fl m (z) — 1 and Q m is the fraction of the critical 
density p c in the form of matter (p m — fl m p c ) [321 ]. 

The initial velocity distribution must be taken to be 
consistent with the mass distribution so that the effects 
of decay are not confused with relaxation of an initially 
nonequilibrium halo to equilibrium. We take the velocity 
ellipsoid to be isotropic, a reasonable first approximation 
to dark-matter halos [33|. This allows us to use Edding- 
ton's formula, 



/(£) 



l 



V8ir 2 



d 2 p 



(4) 



to relate the particle distribution function (DF) f(£) as a 
function of the negative energy £ {£ > for all particles 
bound to the halo) to known functions, the density and 
gravitational potential profile 

We cannot use Eq. ([1} in Eq. Q for arbitrary r; the 
mass of a halo corresponding to a density profile given 
by Eq. (U) is formally infinite. Thus, we must truncate 
the mass profile of the halo, but with care such that the 
distribution function is everywhere positive-definite. We 
use the truncation scheme advocated by Kazantzidis et 
al. (2004) [H, 



p(r) 



r <r r 



M r w (5) 



The density profile has been exponentially truncated at 
the cut off radius r c and with the parameter control- 
ling the sharpness of the transition. Continuity of the 
function and its logarithmic slope demand that 



1 + 3r c /r s 
1 + r c /r s 



I'd 



(6) 



where p s is the scale density, and r s is the scale radius, 



To set the initial positions and velocities of the par- 
ticles, we sample the DF using an acceptance-rejection 
method to find the initial radii and speeds [351 ] - The ra- 
dius and speed are then isotropically mapped to three 
spatial coordinates. 
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A. Parameters 

We begin every simulation with identical initial con- 
ditions, varying only the concentration c of the NFW 
profile, particle lifetimes, and kick velocities. We simu- 
late M v ; r = 10 12 M Q halos only, but we can extrapolate 
our results to other halo masses by considering the struc- 
tural changes to the halos as a function of the decay pa- 
rameters with respect to virial parameters, as discussed 
in Sec. IIII Al To determine the virial spherical top-hat 
overdensity A v , and hence the virial radius, we use a 
cosmology with h = 0.7, Q m = 0.3, and £2a = 0.7. We 
perform a set of 50 simulations, with v-^ — 10, 50, 100, 
200, and 500 km s" 1 and r = 0.1, 1, 10, 50, and 100 
Gyr, and with c = 5 and c = 10. Since the virial speed 
of the halo is u v ir = y/GM V i T / R V i T « 130 km s _1 using a 
spherical top-hat overdensity definition of virial param- 
eters, the kick speeds were chosen to bracket the virial 
value. The typical crossing time of a particle in the halo 
is ~ 500 Myr, and the virial time t vlI = R v i r /v v [ r is a 
few Gyr, meaning that our choices for r bracket the typ- 
ical dynamical time scales of particles in the halo. We 
choose c = 5 and c = 10, since these are typical ACDM 
concentrations for galaxies, groups, and clusters. We use 
observations of such systems to set constraints on the 
decay parameter space in Sec. IIII Bl 

The fiducial simulations used for our analysis in Sec. lIIII 
contain 1.5 x 10 6 particles total, or 10 6 particles within 
the virial radius at the start of the simulations. To test 
the convergence of the simulations, we run the same sim- 
ulations over again but with an order of magnitude fewer 
particles. The convergence tests are further discussed in 

Sec. EES 

We set GADGET-2 parameters, such as the softening, 
maximum time step, and the force accuracy criterion, 
according to the recommendations of Power et al. (2003) 
36]. For the fiducial simulations, the softening length 
is set to 1.3 kpc, and it is set to 3.3 kpc in the low- 
resolution simulations. For most of the simulations, the 
maximum allowed time step is set to Ai max = 10 Myr, 
in accordance with Power et al. (2003). However, for 
simulations with r = 100 Myr and r = 1 Gyr, At max = 1 
Myr such that P max <C 1. We run the simulations for 10 
Gyr. 

We also performed simulations without decay to check 
that our equilibrium initial conditions were indeed in 
equilibrium, and to determine the extent of numerical 
relaxation at the halo centers. We found that the halos 
were stable over the 10 Gyr duration of the simulations, 
and that numerical relaxation ceased to be important for 
r > 4 kpc for the fiducial simulations and r > 10 kpc for 
the low-resolution simulations. 



III. RESULTS 

In Sec. IIII Al we present and summarize the results of 
our simulations. We show how the halo mass, NFW con- 



centration, and NFW goodness of fit depend on the decay 
parameters «k and r, initial concentration, and time. In 
Sec. IIII Bl we show how the results of the simulation can 
be used in combination with the observed galaxy-cluster 
mass function and the halo mass-concentration relation 
to constrain the decay parameter space. 

A. Halo profiles as a function of vt and r 

We present the halo density profiles, multiplied by the 
radius to highlight differences, as a function of Vk, r, and 
time in Fig. [TJ We show the density profiles for halos 
with initial c = 5 since the qualitative features of the 
c = 5 and c = 10 halos are similar, and omit figures for 
t = 0.1 and 1 Gyr with — 200 and 500 km s" 1 since 
the halos are effectively destroyed within a few Gyr. The 
higher lines in each subplot correspond to earlier times, 
and the lines go (top to bottom) t = 0, 2.5, 5, 7.5, and 
10 Gyr after the start of the simulation. 

There are other interesting features of note in Fig. [T] 
First, we see that for small i>k relative to the virial speed, 
there is little change to the dark-matter profile. This 
is expected, as such small «k leads to neither significant 
mass loss nor kinetic energy injection in the halo. How- 
ever, for «k = 50 km s _1 (vk/iVh- ~ 0.4), changes to halo 
structure begin to become significant for r < 10 Gyr. In 
general, changes to halo mass and halo structure are most 
extreme for large and small r. Changes are abrupt 
when r is small, although it takes ~ 5 Gyr for the halos 
to reach equilibrium if r < 1 Gyr. 

We see similar trends in the velocity data, which we 
present in Figs. [5] and [3] In Fig. [21 we show the velocity 
dispersion profiles 

as a function of decay parameters and time, and in Fig. 
[31 we show the velocity anisotropy 

where ay , ere, and <j<p are the velocity dispersions in spher- 
ical coordinates. We set up the initial conditions such 
that a r = <jg = <T0 (such that f3 = 0) and given the equi- 
librium initial conditions, the velocity dispersion is ini- 
tially equivalent to the rms velocities. Even at late times 
in the cases of small t or large i>k, the rms velocities are 
only slightly different from the velocity dispersions. 

As expected, the velocity dispersion profiles <r(r) de- 
crease in normalization as more halo particles decay, since 
the negative heat capacity of a self-gravitating system 
ensures that an injection of kinetic energy due to decays 
results in an overall decrease in the equilibrium kinetic 
energy of the halo (cf. [25|). The most extreme drops in 
the velocity dispersion occur for small r or large for 
which the kinetic energy injected into the halos due to 
decays is large. The velocity dispersion profile quickly 
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FIG. 1: The halo dark-matter density profile times the ra- 
dius, as a function of radius and time. The lines correspond 
to different time slices in the simulation. The solid line rep- 
resents the density at the beginning of the simulations; the 
dashed line is the density 2.5 Gyr after the simulation start; 
the long-dashed line is that after 5 Gyr; the dot-dashed line is 
the density profile after 7.5 Gyr; and the long-dash-dotted line 
represents the density profile after 10 Gyr. Density profiles 
for v k > 200 km s" 1 and r < 1 Gyr are not shown because 
the halos are destroyed within a few Gyr. 



settles to a new equilibrium for r < 1 Gyr. When both 
is large and r is relatively small, the velocity disper- 
sion increases as a function of time at large r. These 
regions are no longer within the virial radius and con- 
tain a number of unbound particles, the latter of which 
in particular drives up the velocity dispersion. 

The velocity anisotropy /3 hovers around zero for most 
times and radii of most of the decay parameters we con- 
sider. This is unsurprising, since the initial conditions 
and the decays have a high degree of spherical symmetry 
in both configuration and velocity space. The exceptions 
in /3 ~ occur, unsurprisingly, for small r and large v^. 
In the case of r < 1 Gyr for — 50 or 100 km s , 
j3 < in the outskirts of the halos at late times, indicat- 
ing that orbits in the outskirts of those halos become tan- 
gentially biased. Radial orbits in the outskirts of galaxies 
are easiest to strip since their binding energy is lower at 
fixed radius than more circular orbits. For large i^and r, 
though, the orbits clearly become radially biased in the 
outskirts of the halo. These halos have not settled into a 
true dynamical equilibrium, and there is a significant net 
outflow of particles. This net outflow drives the radial 
bias in the orbits, and is more significant in the outskirts 
because any particle which decays onto an unbound orbit 
in the interior of the halo must pass through the exterior 
regions. Moreover, the particles originating in the out- 



FIG. 2: The three-dimensional velocity dispersion as a func- 
tion of radius and time, tiled according to Vk and r. The lines 
correspond to the same time slices as in Fig. [1] 



skirts are more vulnerable to ejection, so the daughter 
particles originating in the outskirts will also have a sig- 
nificant radial bias. The interior regions of the halo are 
less vulnerable to decay (unless «k 3> f v ir) an d the dy- 
namical time is short relative to the decay time, both of 
which serve to preserve the initial orbital structure deep 
within the halos. 

However, since the dark- matter mass profile may be in- 
ferred observationally but the dark-matter velocities are 
not, we focus our analysis on the mass distribution in 
halos as a function of decay parameters and time. Per- 
haps the most striking feature of Fig. []] is that the NFW 
halo profile, Eq. [TJ is clearly not a good fit to many of 
the profiles. In the way in which we have chosen to dis- 
play the density profile information, rp as a function of r, 
the profile should be flat for r/r s <C 1 and be monoton- 
ically decreasing for larger r. However, there is a clear 
turnover in rp in some cases, most notably in the r < 1 
Gyr cases, but also at late times for r = 10 or 50 Gyr 
and i>k > 100 km s _1 . 

In order to find the goodness of fit of the NFW profile 
quantitatively, as well as the best-fit NFW parameters, 
we bin the particles radially with respect to the center of 
mass of the system, and use two different likelihood esti- 
mation routines: a simple \ 2 model, which should work 
in the limit of a large number of particles in each bin; and 
a Poisson-based likelihood. We treat the uncertainty in 
the number of particles in each bin as y/~N, where N is the 
observed number of particles in a bin. Instead of using p s 
and r s as free parameters in the NFW fit in Eq. ((T|), we 
use M v ir and c. We use 100 bins between r m ; n (= 4 kpc, 
the minimum radius at which numerical relaxation ceases 
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FIG. 3: The velocity anisotropy f) = 1 — (trf + a^,)/2a^ as a 
function of radius and time, tiled according to Vk and r. The 
lines correspond to the same time slices as in Fig. [T] 



to matter) and r max = 250 (the approximate virial radius 
at t = 0). We bin the particles either linearly or loga- 
rithmically in r, and obtain identical fits, as expected. 
The minimum \ 2 parameters are slightly biased relative 
to those obtained with the Poisson likelihood function, 
but not at a significant level. Thus, we can use statistics 
associated with x 2 to describe the fits, and define x 2 per 
degree of freedom (x 2 /d.o.f.) as our metric for the fit 
quality. 

We summarize the best-fit NFW concentrations and 
the goodness of fit (parametrized by x 2 /d.o.f.) in Fig. |4] 
The concentrations are extrapolated to a Hubble time 
<h = 13.7 Gyr from our simulations, and their values 
correspond to the shading for each grid point represent- 
ing a different (v^,t) pair used for a simulation. The size 
of the grid points corresponds to the quality of the fit. 
If the x 2 /d.o.f. exceeds 1.5 at any time during the sim- 
ulation, the grid point is small; if x 2 /d.o.f.< 1.5 for the 
whole simulation (sampled every 0.5 Gyr), then the point 
is large. In some cases, the fit is somewhat marginal 
(x 2 /d.o.f. ~ 1.5) much of the time, in the dividing re- 
gions of Vk — t parameter space between good and really 
bad NFW fits. The general features between the simula- 
tions with initial c = 5 and c = 10 are similar, though; 
for t < 10 Gyr and v k /v vir > 0.4 (u k > 50 km s _1 ), the 
NFW profile is generally not a good fit to the halo profile 
at some (especially late) times. For lower Vk or high r, the 
NFW density profile provides a reasonable fit throughout 
the lifetime of a halo, even if the halo mass and concen- 
tration change. The concentrations drop rapidly as Vk 
increases, and the final concentration is less sensitive, in 
general, into r than to Vk- The concentration always de- 



creases as a result of decay, which is expected if kinetic 
energy is injected to the halo or in the case of adiabatic 
expansion. 

At large r (specifically, r = 50 and 100 Gyr) and 
«k = 500 km s~\ the changes to the halo structure 
with time are well-described by the adiabatic-expansion 
scheme explored in Ref. [25} ■ If r is much greater than 
the halo dynamical time (tdyn ^5 1 Gyr for particles in an 
M v ij. = 10 12 M Q halo) and the kick speed is high enough 
that daughter particles are cleanly ejected from the halo, 
then changes to the gravitational potential are adiabatic 
and easy to model. If all particles are on circular or- 
bits, the density profile remains constant, but with the 
scale radius rescaled as r s = r s (t = 0)/(l — /) and the 
NFW scale density rescaled as p s = (1 — f) 4 p s (t = 0), 
where / is the fraction of the initial dark-matter par- 
ticles that have decayed. We find that NFW profiles 
are a good fit for Vk — 500 km s _1 (vk/^vir ~ 4) and 
r > 50 Gyr with the adiabatic-expansion parameters, 
but this does not work so well for t = 10 Gyr or for 
Vk = 200 km s _1 (vk/iVir ~ 1-5). Thus, it appears that 
t 3> tdyn and Vk/vyir > 1 are required for the analytic 
adiabatic-expansion description to describe halo evolu- 
tion. 

We reported the goodness of fit and concentration in 
terms of i>k/v V ir instead of Vk because we would like to 
extrapolate the findings of our simulations to other halo 
masses. It is natural to parametrize the kick speeds Vk in 
terms of the virial speed; the virial speed is a measure of 
the binding energy of the halo, and Vk is a measure of the 
kinetic energy injected into the halo. Following Ref. [25| . 
we link the decay time r with the crossing time r/u c irc a t 
the half-mass radius, where u c i rc = [GM(r)/rf 2 , such 
that the typical dynamical time of a particle in a halo is 



t dyn a (GA vPm )- 1 / 2 c- 3 / 4 . 



(9) 



Thus, we can extend our simulation results to halos of 
other masses. 

We summarize the mass remaining after a Hubble time 
in in Fig.O In this plot, we show regions in the decay pa- 
rameter space in terms of Vk/v v i T and r/tdyn on the left, 
and Vk/v v i r and r on the right, in order to highlight both 
the absolute and relative time scales. The Hubble time 
is a fixed time scale, and one can consider r in reference 
to that timescale, or think of tn in terms of dynamical 
times, and then consider the relative time scale r/tdyn- 
The long-dashed line shows M vir (t H )/iVf vlr (i = 0) = 0.9 
for the c = 5 simulations, which has been found via inter- 
polation of our grid of M v i r (tH)/M v i r (t = 0) as a function 
of the decay parameters. This line indicates the transi- 
tion from decays having little effect on the virial mass of 
the halo to having a nontrivial effect. The short-dashed 
line shows M v i r (tn)/M V i r (t = 0) — 0.2, which is around 
the transition between decays causing total destruction of 
the halo and moderate mass loss in the halo. The dotted 
and dot-dashed lines show M v i r (t-n)/M V i T (t = 0) = 0.9 
and M v i r {t}i)/M v i r (t = 0) = 0.2 for the halos which orig- 
inally had c = 10. 
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FIG. 4: Concentration at for an initial c = 5 (left panel) and c — 10 (right panel). The shading indicates the NFW 
concentration at tu, and the size of the points indicates the goodness of the NFW fit. Big points show that the NFW profile 
is a good fit to the density profile throughout the simulation, and small points indicate significant deviations from NFW fits 
during the simulations. 



100 



100 



10 



1 



0.1 



0.01 



M vir (t H )/M i vir > 0.9 



0.01 



0.1 



M vir(W/ M i, v ir < 0-2- 



10 



V v v 



10 



0.1 



0.01 



M vir (t H )/Mi „ ir > 0.9 



0.1 



M vir(tHV M i,vir < °- 2 



10 



FIG. 5: Virial mass remaining at z = as a function of the initial concentration, r, and v^. The long- and short-dashed 
lines show where the virial mass after a Hubble time is 0.9 and 0.2, respectively, times the initial virial mass if the initial 
concentration is c = 5. The dotted and dot-dashed lines show the equivalent 0.9 and 0.2 remaining mass fractions for an initial 
c = fO. Left panel: Remaining mass as a function of r/tdyn and Uk/^vir- Right panel: Remaining mass as a function of r in 
Gyr and «k/^vir- 



When the decays are parametrized in terms of r/tdyn 
and Vk/iMr, the M vix (tn) / M viT (t = 0) = 0.9 and the 
M v i r (tn)/M v i r (t = 0) = 0.2 lines are largely similar be- 
tween the two different initial concentrations. The largest 
discrepancy occurs near Wk/u v ir ~ 0.2 and r < 10 Gyr. 
This effect is real, and is due to the fact that more 
concentrated halos have more particles that are more 
tightly bound to the halos, which makes the particles 
harder to eject. However, for r/tdyn > 10 or for the 
M v i r (tii)/M v i T (t = 0) = 0.2 line, the differences in mass 
retention between the two sets of simulations are small. 
The differences are greatest at small r/td yn and fk/iVir 
because the main particles affected and likely to be kicked 
out of the halo due to such decays are in the outer halo. If 
"fk/wvir ~ 1, particles even relatively tightly bound to the 
halo have a high chance of being ejected, which amelio- 
rates some differences due to differences in concentration. 



In general, though, even though the dynamical times of 
the c — 5 and c = 10 halos are different by a factor of 1.7, 
and even though there are different numbers of dynam- 
ical times in in for the two sets of halos, the remaining 
halo mass at in is not a strong function of concentration 
in r/tdyn ^ Wk/u v i r parameter space. 

The differences are more pronounced in the right-hand 
side plot in Fig.[5l in which the remaining mass is plotted 
as a function of r and t>k/u v ir- The c = 10 halos tend to 
retain more mass at fixed i\,i r and r because the typical 
particle has a higher binding energy than in the c = 5 
case. Since the dynamical time depends on c, which is 
related to the distribution of particle energies within the 
halo, some of the dependence of halo mass loss on con- 
centration is divided out by parameterizing mass loss in 
terms of r/td yn instead of r. 

A more direct comparison would be to examine mass 
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loss and concentration changes as a function of v^/v v i T 
and r/tdyn for a fixed number of dynamical times. When 
we do this, we find that the M v [ r (tn) / M V [ T (t = 0) = 0.2 
lines lie practically on top of each other for c = 5 and 
c = 10, but the M v i r (t H )/M vir (t = 0) = 0.9 line is offset 
by approximately a factor of 2 in T/tdyn between c = 5 
and c = 10. However, what matters in comparisons to 
observations is the effect of decay at a fixed time. Instead 
of mapping results in terms of Uk/fvir and r/idyn, we 
interpolate in terms of i>k/^vir, t, and c. Our simulations 
span the typical ACDM concentrations expected for the 
sizes of halos relevant to the observations, and the virial 
speeds span a factor of 10, a relatively small range. 



In summary, we find that for halos that initially have 
c = 5 or 10 and have the equilibrium NFW halo pro- 
files expected in ACDM simulations, the effects of decay 
in the regime r/idyn < 10 and i> v i r /vk ~ 1 are to dras- 
tically lower the halo concentration in a Hubble time, 
drive substantial mass loss, and drive changes to the 
halo profile such that the NFW profile is not a good 
fit. For r/t dyn > 50 Gyr or v k /v vll < 0.2, the NFW 
profile remains a good fit to the density profile, mass 
loss is not too significant, and the concentration does 
not drop too much. Mass loss and concentration reduc- 
tion are more severe for fixed r/idyn for large until 
Vk/vvir ~ 4, at which point the decay can be described 
analytically in the adiabatic-expansion approximation, 
and the effects of decay no longer depend on i>k as long 
as Wk/wvir ^4. At such high i>k, any daughter parti- 
cle in a decay is ejected from the halo, so there is no 
kinetic energy injection into the halo as a result of de- 
cays. In the intermediate regime [roughly corresponding 
to the space between the M v i T (tB.) /M v - lt (t = 0) = 0.2 and 
M v i T (tH)/M v i r (t = 0) = 0.9 lines], mass loss is moder- 
ate although the concentration may drop quite a lot, and 
the inner part of the density profile is shallower than for 
NFW. 



There are some differences between the c = 5 and 
c = 10 simulations beyond that which can be factored 
out by parameterizing decay in terms of the virial time 
scales Wk/wvir and r/idyn- This is due to the fact that the 
gravitational potential, and hence binding energy, does 
depend on the concentration in a way that is not com- 
pletely factored out by the virial scaling. Even though we 
parametrize halos in terms of "typical" dynamical times 
and typical speeds, there is a diversity of particle speeds 
and time scales within a halo. Nonetheless, since our 
simulations span the range of idyn and c of the halo ob- 
servations below, we are able to interpolate our simula- 
tions and map ACDM dark-matter halos to halos in a 
decaying-dark-matter cosmology. 



B. Relation to observational constraints 

1. The mass- concentration relation 

Under the hierarchical picture of structure formation 
in the Universe, it is expected that low-mass dark-matter 
halos form earlier than high-mass halos. As a conse- 
quence of this early formation time, the inner parts of 
low-mass halos are far more dense relative to the out- 
skirts of the halo relative to high-mass halos [37|. In 
other words, low-mass halos are more concentrated than 
high-mass halos. This trend has been observed and quan- 
tified in simulations of ACDM cosmologies [38- 43]. It is 
found that cosmologies with higher matter fractions, f2 m , 
and a higher amplitude of matter fluctuations averaged 
in 8 h~ 1 Mpc spheres, as, produce higher concentrations 
for fixed halo mass than cosmologies in which f2 m and a 8 
are low. This is due to the fact that structures collapse 
earlier when the Universe is more dense if Q m and a& are 
higher. 

Decays will make halos appear less concentrated; 
hence, f2 m and as inferred from halos that have expe- 
rienced decay if one begins with a WIMP dark-matter 
ansatz would be lower than that inferred from observa- 
tions of earlier epochs. We can set conservative con- 
straints on the decay parameter space by comparing ob- 
servationally inferred mass-concentration relations with 
the mass-concentration relation derived from decays in a 
high-f2 m , high-erg cosmology. The CMB last-scattering 
surface is the earliest window into cosmological parame- 
ters. The six-parameter cosmological fit to the Wilkinson 
Microwave Anisotropy Probe (WMAP) seven-year data 
set suggests fl m = 0.267 ± 0.026 (1-a) and a 8 = 0.801 ± 
0.030 (1-a) 0. The 2-a upper limits of Q m and a 8 
correspond approximately to the mean WMAP first-year 
mean parameters (45| . The mean mass-concentration re- 
lation found in simulations of a ACDM cosmology using 
the WMAP-1 parameters is illustrated in Fig. [51 and is 
denoted by the upper solid line [42| . 

We use our grid of mass loss and concentration as a 
function of r and i>k/v V ir, found using the simulations 
and described in Sec. IIII A[ to map ACDM halo masses 
and concentrations in a WMAP-1 cosmology to those of 
halos at t = %. We show the mass-concentration relation 
as a function of decay parameters in Fig. [6j along with 
the mean theoretical and observed relations. It should 
be noted that we use the best-fit NFW concentrations 
[Eq. ©] in Fig. |H1 even if the NFW profile is not a good 
fit to the density profile. Nonetheless, the NFW con- 
centration provides a reasonable measure of the density 
contrast within the halo. 

In Fig. |6l the solid lines represent the mean mass- 
concentration relations found in dissipationlcss, baryon- 
free cosmological simulations using WMAP one-year 
(upper line) [45[ and WMAP three-year (bottom line) 
[46j | cosmologies (42J. The long-dashed lines represent 
decaying-dark-matter models with (top to bottom) t — 
100, 40, 20, and 5 Gyr under the assumption that the 
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WMAP one-year cosmology is the underlying cosmology. 
The short-dashed and dotted lines are observationally in- 
ferred mean mass-concentration relations using a sam- 
ple of clusters observed in a variety of ways [47[ and x- 
ray observations [48|, respectively. The error bar in the 
upper left corner of the plots shows the scatter in the 
mass-concentration relation for simulations, and the typ- 
ical uncertainty in the mean mass-concentration relation 
from the observations. The shaded region corresponds to 
the 1-a region of the mean mass-concentration relation 
inferred from weak-lensing observations of Sloan Digital 
Sky Survey (SDSS) galaxies [U. 

We use only the weak-lensing mass-concentration rela- 
tion to constrain the decay parameters. This is because 
the cluster and x-ray samples have serious selection ef- 
fects which have not been fully quantified. Reference 47 1 
finds that the concentration of strong- lens clusters is even 
higher than that predicted by Ref. [51], who identified 
strong-lens systems in TV-body simulations. Dark-matter 
profiles can only be inferred from x-ray observations if the 
halos are relaxed; relaxed halos have, on average, higher 
concentration than the population of halos as a whole 
39]. Thus, until the selection effects are better quanti- 
fied for the x-ray and strong-lens systems, one should not 
use the mass concentration inferred from observations of 
those systems to constrain dark-matter models. 

Even though there are also systematics that have not 
been fully quantified in the weak-lensing measurement 
(see Sec. HID of Ref. [49j] for a sampling of possi- 
ble effects), the measurement provides conservative con- 
straints because it lies slightly below simulation values, 
on which we base our decay predictions. We consider 
a (vk,r) pair of parameters to be excluded if its mass- 
concentration relation lies below the 1-a lower limit of 
the mean weak-lensing mass-concentration relation. 

Using this criterion, we see from Fig. [5] that current 
weak-lensing measurements are unable to constrain r for 
Vk ^ 100 km s _1 . However, r < 30 Gyr is ruled out for 
300 km s _1 < v k < 700 km s _1 , and r < 35 - 40 km s" 1 
for larger values of v^. Interestingly, the constraints from 
the simulations for large are nearly identical to those 
found analytically using a simple model of adiabatic ex- 
pansion (valid only if r/tdyn 3> 1 and Uk/iW ^ 1) in 
Ref. [H . 

We have neglected the effects of baryons on the halo 
concentration. This is largely because it is not clear 
how baryons affect the dark-matter halo. In cosmo- 
logical simulations with two different prescriptions for 
baryon physics, Rudd et al. (2008) [52| found either lit- 
tle difference relative to dark-matter-only simulations, or 
an increase in concentration of ~ 50% above the dark- 
matter-only values. However, even in the latter case, 
the minimum allowed r for fixed i>k would only decrease 
~ 15% - 30%. 

Better constraints may be obtained in the future. Up- 
coming deep wide-field surveys will increase the statistics 
of weak and strong lensing and x-ray cluster profiles, and 
probe down to smaller halo masses [53Tj56| . However, 



improvements in using the (redshift-dependent) mass- 
concentration relation to constrain dark-matter decays 
(or cosmological parameters) will require a better un- 
derstanding of the systematics and selection effects for 
the various methods of measuring the dark-matter halo 
masses of galaxies and clusters. This is likely to be a sub- 
ject of significant future work, since the selection effects 
are relevant to determining halo mass functions. The 
evolution of halo mass functions depends on the growth 
function. Since the growth function is sensitive to the 
equation of state of dark energy, there will be signifi- 
cant efforts to understand any systematic errors that may 
obscure inferences about dark energy from observations 
[57I . In addition, there needs to be a better theoreti- 
cal understanding of which processes of galaxy evolution 
drive changes to dark-matter halo profiles. This sub- 
ject has generated recent interest [58l462j |. and with con- 
tinuing improvements in the spatial resolution and star- 
formation and feedback recipes in simulations, it should 
be possible to determine which physical processes influ- 
ence the distribution of dark matter within galaxies. 



2. Galaxy-cluster mass functions 

The mass function of galaxy clusters n(> M v i r ) is sen- 
sitive to il m and as for similar reasons as for the mass- 
concentration relation; if the amplitude of matter fluctu- 
ations is higher, structure forms earlier and is more abun- 
dant in the local Universe [6314671 ]. Since galaxy clusters 
are the largest and rarest of virialized structures, corre- 
sponding to small fluctuations in the linear matter den- 
sity field smoothed on ~ 10 Mpc scales, small changes 
in as and f2 m will result in large changes to the cluster 
mass function. 

If dark matter were to decay into relativistic particles, 
the cluster mass function would look different than the 
ACDM mass function because the continuous pumping 
of relativistic particles into the Universe would alter the 
background evolution of the Universe and hence the den- 
sity threshold for halo collapse. This effect was inves- 
tigated by Oguri et al. (2003) [68|. However, the main 
effect of decays in which the mass splitting between dark- 
matter species is small is to reduce the mass of halos due 
to ejection resulting from decays. Since as and Q m should 
be the same regardless of the redshift of the matter tracer 
in a ACDM universe, the smoking gun for decays would 
be if fi m and as inferred from an early epoch (the cosmic 
microwave background or Lyman-a forest) were system- 
atically higher than that inferred from the z — cluster 
mass function. 

We use parameter estimates from the WMAP seven- 
year data set and low-redshift observations of clusters to 
constrain decay parameter space. As with the analysis 
using the mass-concentration relation, we find conserva- 
tive bounds on the decay parameter space by considering 
the parameters for which the cluster mass function with 
f2 m and as set to 2-er above the WMAP-7 parameters 
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FIG. 6: The mass-concentration relations predicted for ACDM and decaying-dark-matter cosmologies and inferred from obser- 
vations of galaxies, galaxy groups, and clusters. The short-dashed line shows the mean mass-concentration relation found in 
Ref. [50(| using cluster observations, and the dotted line shows that from x-ray observations in Ref. [481 ]. The shaded region shows 
the To" fits to the mean mass-concentration relation inferred from weak lensing [4S| . The upper and lower solid lines show the 
mean mass-concentration relation found in dissipationless simulations using WMAP-1 and WMAP-3 cosmological parameters, 
which bracket the WMAP-7 mass-concentration relation [42|. The long-dashed lines show the mass-concentration relation for 
decaying dark matter, with the lines corresponding to (from thickest to thinnest) r = 5, 20, 40, and 100 Gyr. The error bar 
shows the intrinsic scatter in the theoretical relations, and is also approximately the 1-cr error in the mean observationally 
inferred relations. 



maps to a cluster mass function that is 1-cr below f2 m 
and us inferred from observations of low-redshift clusters. 
Both optical and x-ray observations of clusters indicate 
that the 1-a lower limits on fi m and as are approximately 
the same as the 2-er lower limits on those parameters from 
WMAP-7, fl m w 0.2 and a 8 « 0.7 d i, MEM ■ 

We use the relation between the ACDM halo masses 
and the halo masses after t = in with decay to find the 
cluster mass function, 



n f (> M f ) 



dMi 



,dni(Mi(M f ')) dMi 



Mf 



dM-, 



dMi 



(10) 



Here, rif (> Mf) is the comoving number density of ha- 
los with masses larger than Mf after a Hubble time of 
decay, dn-JdM-, is the Tinker et al. (2008) [73[ ACDM 
differential mass function, M; is the initial ACDM mass 
of the halos, and Mf is the final halo mass after a Hub- 
ble time of decay. Since the relation between M; and Mi 



masses depends on concentration (Sec. MI Al . we use the 
WMAP-1 mass-concentration relation from Ref. [42[ to 
set the initial concentrations of the initial ACDM halos 
and interpolate among our simulations to find Mf from 
M. 

We show the galaxy-cluster mass function as a function 
of decay parameters in Fig. [7] The dotted line is the mass 
function for the high-f2 m , high-os cosmology, the solid 
line is for the mean WMAP-7 cosmology, and the long- 
dashed line represents the mass function for a cosmology 
set to the 2-er lower limits of Q m and as parameters from 
WMAP-7 from the six-parameter cosmological fit [iij . 
All other cosmological parameters have been set to the 
WMAP-7 mean values, as those parameters do not affect 
the normalization or shape of the mass function nearly 
as much as f2 m and as- The short-dashed lines show the 
mass functions assuming (top to bottom) r = 100, 40, 
20, and 5 Gyr. Each panel shows a different i>k- 

It is clear that the cluster mass function (for which 
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FIG. 7: z=0 cluster mass functions. In each panel, the dotted (blue) line represents the Tinker et al. (2008) [z3] mass function 
for a cosmology with Q m = 0.318 and as = 0.868, but with all other parameters set to the mean values found in the six- 
parameter fit of the WMAP-7 data The solid (green) line represents the WMAP-7 six-parameter mean cosmology, and 
the long-dashed (red) line represents a cosmology with fl m = 0.198 and as = 0.724. The short-dashed (black) lines show 
cluster mass functions for a cosmology with Q m — 0.318 and as = 0.868, but with decaying dark matter with (bottom to top) 
r = 5, 20, 40, and 100 Gyr. For clarity, each panel shows a different tv 



we consider M vir > lO 14 Af to be classified as a clus- 
ter) is relatively insensitive to small kick speeds, «k ^ 
500 km s _1 . This is because the typical virial speed of 
a cluster is v v - lr > 600 km s _1 . From Fig. [51 it is ap- 
parent that severe mass loss in dark-matter halos re- 
quires Wk/fvir ^ 1 and t < 10 Gyr. However, the 
shape of the mass function is somewhat different than 
those predicted in ACDM cosmologies, since higher mass 
halos are less affected by lower t>k than low-mass ha- 
los. For «k = 1000 km s _1 , the number density of 
M V j r > 10 14 M Q is approximately the same for r as the 
lowest number density allowed by low-redshift observa- 
tions of clusters. However, even for higher r, the cluster 
mass function is far less steep than predicted by ACDM. 
At ^k = 2000 km s _1 , the cluster number density is too 



low for r < 30 Gyr. At even higher v^, the cluster mass 
functions converge to the analytic estimates of Ref. 25 1 
in the case of Uk/iVir ^ lj an d constrain r > 35 — 40 
Gyr. The shape of the cluster mass function as well as 
the number density n(> 10 14 M Q ) can be used to con- 
strain decays, but the statistics are generally lower for 
higher cluster masses. 

Constraints on dark-matter decays from the cluster 
mass function are more robust to baryon physics than 
the mass-concentration relation. Using iV-body hydrody- 
namic simulations with two different models for baryon 
physics, Rudd et al. (2008) [52[ found that the normal- 
ization of the halo mass function deviates by 10% with 
the inclusion of baryons relative to the halo mass function 
found in A^-body cold-dark-matter-only simulations. 
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FIG. 8: Summary of the mass-concentration and z — cluster 
mass function limits on the decay parameter space. 



The redshift evolution of the cluster mass function is 
also a useful probe of decays [25j]. The signature of de- 
cays will be that, for a fixed z = mass function, the 
high-redshift mass functions will have higher normaliza- 
tion and a different shape than in a ACDM universe. If 
one were to estimate fi m and as for different redshift 
slices with the ansatz of a ACDM universe, one would 
estimate systematically higher as and f2 m at higher red- 
shifts. Currently, as and f2 m inferred from small samples 
of x-ray clusters up to z = 0.7 are consistent at differ- 
ent redshift slices, although the error bars are still fairly 
large [1, |(5{| [Zll uM ■ Upcoming dee p optic al, x-ray, and 
Sunyaev-Zel'dovich surveys |53T - l5q , |74h76I ] should allow 
for better constraints on decays, since they will find large 
samples of clusters at a variety of redshifts. 



C. Summary 

We summarize the constraints on the decaying-dark- 
matter parameter space in Fig. [8l The exclusion region, 
marked "ruled out" on the plot, comes from three types 
of observations. For relativistic v^, the WMAP-1 data set 
the strongest constraint on r, with r > 123 Gyr at 68% 
C.L. [77|. This restriction is essentially a restriction on 
p r (z), the radiation density of the Universe as a function 
of time. For v k }t 2000 km s _1 , both the galaxy-cluster 
mass function and the mass-concentration relation re- 
quire r > 40 Gyr. At lower Vk, the mass-concentration 
relation sets best constraints on r, at approximately the 
level r > 30 Gyr for 300 km s _1 < v k < 700 km s" 1 and 



t > 20 Gyr for 100 km s 



< 300 km s" 



We note that this exclusion region looks similar to that 
of Ref. [25] . However, the exclusion region in this paper is 
more robust; the exclusion region for Vk < 5000 km s _1 
in Ref. [25| was based on noting that the clustering of 
galaxies in the SDSS appeared to be consistent with the 



clustering of ACDM halos with M vir > 10 i2 M Q [78j. 
The constraint was set based on not letting ACDM halos 
with M v i r ~ 10 12 Mq lose more than half their mass due 
to decays. This is a fairly hand- waving constraint. It 
will likely take cosmological simulations of decaying dark 
matter in order to map galaxy clustering to constraints 
on dark-matter decays. 

The regions labeled "CDM-like" ( "cold-dark- matter- 
like", such as WIMPs) and "allowed" are still plausible 
regions of parameter space. However, it will be chal- 
lenging to probe the region marked "CDM-like". For 
t > 100 Gyr, there is hardly any change to halo proper- 
ties; mass loss relative to ACDM halos is on the ~ 10% 
level, and concentrations go down by at most ~ 20%. 
The NFW profile is still a reasonable fit to the halo den- 
sity profiles. It is not clear if the systematics in even the 
largest and deepest all-sky surveys will be small enough 
to probe such subtle deviations from ACDM halo prop- 
erties. Decay models with Vk < 1 km s _1 will also be 
challenging to probe, as such speeds are far less than the 
typical virial speeds of galaxies or larger systems. The 
smallest halos that have been found observationally have 
M vir ~ 1O 9 M , and are subhalos at that [79j | 80t . S uch 
halos have iVir ~ 13 km s _1 . As we saw in Sec. lIII Al the 
structural properties of dark-matter halos hardly change 
if Vk/v vir <C 1. Moreover, there are poor statistics of 
such small halos, and it is not clear how many halos of 
that size will host a luminous component. The statistics 
of such small objects may be probed using gravitational 
milli-lensing, but we are years away from building any 
sort of statistical sample of such halos [8ll - l83l |. 

There is some hope of further constraining the re- 
gion marked "allowed" in the near future. Sunyaev- 
Zel'dovich, x-ray, and deep optical surveys should un- 
cover ~ 10 5 galaxy clusters up to quite high redshift [53r - 
HEIzEIZSl]; then, both the z = cluster mass function and 
the evolution of the cluster mass function could be used 
to more accurately constrain the decay parameter space. 
The deep optical surveys should also allow for a more 
precise measurement of the mean mass-concentration re- 
lation, and down to lower halo masses. On the theoretical 
side, a better understanding of selection effects will allow 
for a better comparison of the decaying-dark-matter pre- 
dictions with observations. In addition, if cosmological 
simulations of decaying dark matter are performed, one 
can hope to characterize the halo occupation distribution 
(HOD) [78l. [84T - l88l |. which will open up galaxy clustering 
as an avenue to better explore the decay parameter space. 

So far, we have not explored using dwarf galaxies to 
constrain the allowed region of parameter space, as sug- 
gested by Ref. 24]. This is because there are no robust 
measurements of the abundance of M vu < IO^Mq ha- 
los, nor are there strong constraints on the inner pro- 
files of dwarf galaxies [89-92]. 1 One could imagine that 



1 We also disagree with several of the technical aspects of Ref. |24| : 
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the M v i r ~ 10 9 M Q dwarf galaxies could be the result 
of decays of halos born with much higher virial masses, 
although the density profile would be significantly differ- 
ent than expected in ACDM. In the future, large optical 
surveys such as LSST will look farther down the lumi- 
nosity function, both in the Local Group and in the local 
universe [55l |93| , which will allow for better constraints 
using the abundance of low-mass halos. Astrometric ob- 
servations of stars in Local Group dwarf galaxies will be 
critical getting the three-dimensional stellar kinematic 
velocity data for better constraints on the dark-matter 
density profiles in those galaxies [90]. This will be pos- 
sible with next-generation astrometric satellites, such as 
Gam and SIM MM. 



IV. DISCUSSION 

In the previous section, we presented the results of our 
simulations and showed how those results, along with 
the mass-concentration relation and cluster mass func- 
tions inferred from existing observations, could be used 
to constrain the decaying-dark-matter parameter space. 
In this section, we address two issues that will affect fu- 
ture attempts to further constrain the decay parameter 
space: numerical convergence in simulations, and more 
accurate mapping of simulations to observations. 



A. Resolution effects 

We based all our analysis in Sec. [TTT] on simulations 
that initially had ~ 10 6 particles within the virial radius. 
We would like to know which of the results we described 
above are fully converged, and which results may be spu- 
rious due to resolution effects. Moreover, we would like 
to understand resolution effects before embarking on cos- 
mological simulations, which are far more costly than the 
simulations described in Sec. [IT] If > 10 6 particles are re- 
quired within a virial radius in order to determine halo 
properties, it becomes prohibitive to get good statistics 
on dark-matter halos (and subhalos) on a variety of mass 
scales in cosmological simulations. 

In order to understand the robustness of our results, 
and to determine if one can get away with fewer parti- 
cles in the virial radius for cosmological simulations, we 
perform simulations with the same c, t>k, and r values as 
those used in Sec. IIII1 but with 10 5 particles within the 
virial radius. The low-resolution simulations are briefly 
described in Sec. [XT] Our goal is to determine how the 
virial mass, concentration, and goodness of fit to the 



(1) The mechanical kinetic energy injected by the kick velocity 
is far more important than the change in the rest-mass energy 
they consider. (2) Our analysis shows that the preservation of 
a NFW profile, assumed therein, is not valid. (3) We find that 
halos may undergo disruption even for < v v - lT . 
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FIG. 9: Density profile for the c = 10, v k = 100 km s"\ 
t = 10 Gyr simulations at 9 Gyr divided by the best-fit NFW 
density profile. Solid line: 10 6 particles in the virial radius. 
Dashed line: 10 particles within the virial radius. 



NFW density profile compare between the two sets of 
simulations. 

We find that the virial mass derived from fits to the 
NFW profile are consistent between the two sets of sim- 
ulations, and are consistent with the virial mass as de- 
termined by finding the average density of dark-matter 
particles within a sphere centered on the halo center of 
mass. Similarly, the NFW concentration parameter from 
the fits is also consistent between the two sets of simula- 
tions. Therefore, if one is concerned with accurately de- 
termining halo masses or NFW concentration parameters 
in a simulation with a decaying-dark-matter cosmology, 
one requires no more than 10 5 particles within the virial 
radius. 

By contrast, the quality of the NFW fit can be sensitive 
to the number of particles within the virial radius. For 
example, x 2 /d.o.f. < 1.4 for fits to all particles within the 
virial radius for = 100 km s _1 and r = 10 Gyr for 10 5 
particles within the virial radius, regardless of the initial 
concentration. However, the x 2 is bad for t > 0.5 Gyr for 
both c = 5 and c = 10 if there are initially 10 6 particles 
in the virial radius, even though the NFW fit parameters 
of the two different-resolution simulations are identical 
(see Fig. [TU] for x 2 m the c = 10 case). This is illustrated 
in Fig. [51 in which we show the density profile from the 
simulation (p num ) divided by the best-fit NFW profile 
for the two simulations at t = 9 Gyr with c = 10, «k 
I i Hi km s - , and t — 10 Gyr (pat)- The innermost points 
(r < 10 kpc) show evidence of numerical relaxation. For 
the rest of the halo, one can tell that there is significant 
curvature in p n um/pfit for the 10 6 -particle simulation, but 
not for the 10 5 -particle simulations. This difference is 
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likely in large part due to the increased Poisson noise in 
the 10 5 -particle simulation. In general, the discrepancy 
in the quality of the fit appears to be exacerbated in cases 
in which r < 10 Gyr and Vk/«vir ^ 1- 

The consequences for our results are as follows. Our 
constraints on the decay parameter space from cluster 
mass functions depended on the mapping between the 
ACDM halo mass and the halo mass after a time £h- The 
mass estimated in the high-resolution simulations used in 
Sec. Mil is virtually identical to that found in the lower- 
resolution simulations, so we believe the halo mass (as a 
function of decay parameters and time) to be converged 
and reliable. Similarly, the use of the mass-concentration 
relation depends on the mapping of ACDM halo masses 
and NFW concentrations to those of halos that have ex- 
perienced decays. Again, we believe the virial masses and 
NFW concentrations to be converged in our simulations. 
The shape of the dark-matter density profile is not used 
for mapping either of the observations to constraints on 
decay parameter space. Thus, we believe the mapping of 
our simulations to observational constraints to be robust. 

The implications for future cosmological A-body sim- 
ulations are as follows. We believe it is possible to get ac- 
curate estimates of halo masses and NFW concentrations 
using no more than ~ 10 5 particles in the virial radius. It 
may be possible to get good statistics using fewer parti- 
cles, but it will likely take convergence tests in fully cos- 
mological A^-body simulations to determine the minimum 
particle number required for convergence, especially for 
subhalos. This means that it may be possible to explore 
decaying-dark-matter cosmologies with reasonable statis- 
tics using modest computational resources. The only sit- 
uation in which higher particle numbers may be required 
is a detailed comparison of the shape of decaying-dark- 
matter halo profiles with respect to ACDM halos. In this 
case, we showed that ~ 10 6 particles are required within 
the virial radius to determine the goodness of the NFW 
fit, at least down to ~ 0.02i? vir for Vk/iVir an( l r ^ 10 
Gyr. In other regions of the decay parameter space, the 
10 5 - and 10 6 -particle simulations converge on the good- 
ness of the NFW fit. If one is mainly concerned with 
radical differences between ACDM and decaying-dark- 
matter halo density profiles, then ~ 10 5 particles within 
the virial (or tidal truncation) radius is likely sufficient; 
more subtle differences require higher particle numbers. 



B. Mapping simulations to observations 

In Sec. IIII1 we used NFW fits to the entire halos to re- 
late halo properties to observables. However, astronom- 
ical probes of the gravitational potential of galaxies and 
their halos are rarely sensitive to the entire halo. Rota- 
tion curves rarely extend beyond the inner 10% — 20% of 
the virial radius of the halo 95 97] , strong lensing is gen- 
erally also most sensitiv e to t he mass within ~ 10% — 20% 
of the virial radius [98l-[l00j , and weak lensing is gener- 
ally most sensitive to the outer parts of halos (excluding 



the inner ~ 20% of the virial radius) [9^, Il0l| . X-ray gas 
profiles in galaxies, groups, and clusters can ex tend up t o 
several tens of percents from the halo centers [102L 1 1 031 ] . 
If dark-matter halos were well described by the NFW pro- 
file and there were no systematics in the observations that 
would bias the determination of the dark-matter profile, 
then all probes should yield the same profile fits. 

Different probes sometimes yield different profile fits, 
though. For example, it has been claimed that NFW fits 
to dwarf-galaxy rotation curves trend lower in concentra- 
tions than predicted in dissipationless cold-dark-matter 
simulations, and that NFW profiles are not always good 
fits to the data [42|, [93] . Strong-lens systems and x- 
ray clusters appear to have higher concentrations than 
the average concentrations predicted by simulations, al- 
though a larg e part of that may be due to selection ef- 
fects [HTl Il04j . In some cases, multiple probes are used 
on the same system and show different fits. For exam- 
ple, Ref. [9^] found moderate tension between the weak- 
lensing and strong- lensing NFW fits to Abell 611 (al- 
though the mass contained within 100 kpc, which is ap- 
proximately the transition point between the domains of 
the probes, is consistent), and strong tension with NFW 
fits derived from stellar kinematic data, so much so that 
density profiles significantly shallower than NFW are pre- 
ferred if the three data sets are combined [99l ]. 

We would like to explore the possibility that these dis- 
crepancies could be explained in the context of decaying 
dark matter. In Sec. IIII Bi we showed that the NFW 
profile is not always a good fit to dark-matter halos af- 
ter some of the X particles have decayed, especially if 
fk/^vir ~ 1- Thus, profile fits to subsections of the halo 
may deviate from fits to the whole halo. 

To demonstrate trends in the fit parameters as a func- 
tion of where in the halo one performs the fit, we show 
in Fig. [TO] NFW fit parameters and the goodness of fit 
(x 2 /d.o.f.) as a function of time for a halo initially 
in equilibrium with M v - lr = 10 12 M@ and c — 10, with 
t>k = 100 km s _1 and r — 10 Gyr. The solid (blue) line 
shows the NFW fit parameters M V - 1T and c using all par- 
ticles within the virial radius. The error bars indicate the 
90% range of parameter fits based on bootstrap resam- 
pling particles at each time slice. The short-dashed line 
shows the best-fit parameters using only those particles 
in the inner half of the virial radius, and the long-dashed 
line shows the fits to particles within the inner 10% of 
the virial radius. In these three cases, the inner cut-off 
to the fits is 4 kpc, corresponding to the minimum radius 
at which numerical effects on the density profile can be 
ignored. The dot-dashed line shows the fit resulting from 
considering only those particles with radii > 0.5i? V i r - In 
many cases, the NFW profile does not provide a good fit 
to the density profiles, but it is also not always a good fit 
to observational data. The trends we discuss below are 
there regardless of the goodness of the NFW fit. 

We find that the concentration derived from the inner 
parts of the halo is systematically lower than the concen- 
tration derived from the halo as a whole at late times, 
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FIG. 10: NFW fits to a simulation of a 1O 12 M , c = 10 halo 
with Uk = 100 km s _1 and r = 10 Gyr. The solid line repre- 
sents a fit using all data down to the inner cut-off (4 kpc), the 
short-dashed line represents fits using particles within half of 
the initial virial radius, the long-dashed line shows fits within 
the inner tenth of the initial virial radius, and the dot-dashed 
line shows fits to the outer half of the halo. The error bars 
are the 90% limits from bootstrap resampling particles at each 
timestep, although the error bars lose meaning as the \ 2 be- 
comes large. Top panel: NFW fits to the virial mass. The 
true virial mass follows the solid line. Middle panel: The 
NFW concentration. Bottom panel: x 2 /d-0.f. for the NFW 
fit. 



although at early times, the inner concentration is some- 
times higher for r <C idyn or r ^> t^. However, if the 
maximum radius used in the fit is ~ 0.1i? V j r , the concen- 
tration from the NFW fit is always lower than the con- 
centration from the NFW fit to the halo as a whole. In 
general, the discrepancy among concentrations becomes 
less severe as a larger fraction of the halo particles are 
used for the fit. The NFW fit becomes better and con- 
centrations tend to be slightly higher if the minimum 
radius is raised to 10 kpc from 4 kpc, and in general, 
NFW fits become better as the maximum radius consid- 
ered in the fit decreases. Fits to the outer part of the 
halo show a trend that the halo concentration is initially 
higher than that estimated from the halo as a whole, then 
dips low, and then becomes relatively high again. This 
phenomenon is due to the finite time it takes for parti- 
cles ejected from the inner halo to reach the outer halo, 
and the longer dynamical (and hence, equilibration) time 
in the outer halo. The NFW profile is almost always a 
reasonable fit to the data with r > 0.5i? v ir- 

Are these trends consistent with observations? The 
relatively low concentration measured in the innermost 
part of the halo relative to the halo as a whole is qual- 



itatively consistent with the low concentrations inferred 
from low-surface-brightness and dwarf-galaxy rotation 
curves. However, in the absence of selection biases, one 
would expect the weak-lensing-derived concentrations to 
be slightly higher than strong-lensing-derived concentra- 
tions for the same halos, and also slightly higher than 
concentrations inferred from x-ray gas profiles, which is 
the opposite of what is observed. X-ray and strong- 
lensing halo samples are likely highly biased towards 
higher concentrations, and so a better understanding of 
selection effects is necessary in order to understand the 
differences in inferred concentration among the various 
lcnsing, x-ray, and dynamical probes of dark-matter den- 
sity profiles. 

Nonetheless, the prediction from decaying-dark-mattcr 
models is that it is possible that different probes of dark- 
matter halo profiles will indicate different concentrations 
for the same halo, even if the NFW profile is shown to be 
a good fit to data from probes sensitive to the innermost 
or outermost parts of the halo, such as weak lensing or 
rotation curves. 



CONCLUSION 



In this work, we have used simulations of two-body 
dark-matter decay in isolated halos initially in equilib- 
rium to characterize the response of halos as a function 
of the center-of-mass kick speed Vk, the decay time r, 
and time. We find that for the union of the parameter 
space Wk/iVir 0.2 and r/td yn > 100, dark-matter ha- 
los are essentially unchanged; the mass loss is < 10% 
over a Hubble time, and the shape of the density pro- 
file does not change, although the concentration declines 
somewhat for v^/v v i T > 1. For the union of v\,/v v - lr > 1 
and r/tdyn ^5 10, there is severe mass loss and a radi- 
cal change to the shape of the dark-matter density pro- 
file. For the rest of the decaying-dark-matter parameter 
space, there is moderate mass loss and significant de- 
viations of the shape of the dark-matter density profile 
relative to the initial halo properties. 

We used our simulations in conjunction with the 
observed mass-concentration relation and the galaxy- 
cluster mass function to constrain the two-parameter 
decay parameter space. We find that the mass- 
concentration relation yields the stronger constraint on 
t for a broad range of , constraining r > 30 — 40 
Gyr for v k > 200 km s _1 . Constraints on the decay 
parameters should improve greatly in the next 5-10 
years with blossoming wide-field Sunyaev-Zcl'dovich and 
optical surveys, and with cosmological simulations of 
decaying dark matter. 



15 



Acknowledgments FG03-92-ER40701 and the Gordon and Betty Moore 

Foundation. 

We thank Andrew Benson for helpful discussions. This 
work was supported at Caltech by DOE Grant No. DE- 



[1] R. Kessler, A. C. Becker, D. Cinabro, J. Vanderplas, 
J. A. Frieman, J. Marriner, T. M. Davis, B. Dilday, 
J. Holtzman, S. W. Jha, et al., Astrophys. J. Suppl. 
Ser. 185, 32 (2009). 

[2] B. A. Reid et al., arXiv:0907.1659. 

[3] A. Vikhlinin, A. V. Kravtsov, R. A. Burenin, H. Ebel- 
ing, W. R. Forman, A. Hornstrup, C. Jones, S. S. Mur- 
ray, D. Nagai, H. Quintana, et al., Astrophys. J. 692, 
1060 (2009). 

[4] E. Komatsu, K. M. Smith, J. Dunkley, C. L. Bennett, 
B. Gold, G. Hinshaw, N. Jarosik, D. Larson, M. R. 
Nolta, L. Page, et al., arXiv:1001.4538. 

[5] E. Rozo, R. H. Wechsler, E. S. Rykoff, J. T. Annis, M. R. 
Becker, A. E. Evrard, J. A. Frieman, S. M. Hansen, 
J. Hao, D. E. Johnston, et al., Astrophys. J. 708, 645 
(2010). 

[6] G. Jungman, M. Kamionkowski, and K. Griest, Phys. 

Rep. 267, 195 (1996). 
[7] T. Appelquist, H.-C. Cheng, and B. A. Dobrescu, Phys. 

Rev. D 64, 035002 (2001). 
[8] H.-C. Cheng, J. L. Feng, and K. T. Matchev, Phys. Rev. 

Lett. 89, 211301 (2002). 
[9] G. Servant and T. M. P. Tait, New J. Phys. 4, 99 (2002). 
[10] J. Hubisz and P. Meade, Phys. Rev. D 71, 035016 

(2005). 

[11] B. Moore, S. Ghigna, F. Governato, G. Lake, T. Quinn, 
J. Stadel, and P. Tozzi, Astrophys. J. 524, L19 (1999). 

[12] S. Davidson, S. Hannestad, and G. Raffelt, J. High En- 
ergy Phys. 05, 003 (2000). 

[13] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, 
and M. Kamionkowski, Phys. Rev. D70, 083501 (2004). 

[14] D. N. Spergel and P. J. Steinhardt, Phys. Rev. Lett. 84, 
3760 (2000). 

[15] S. S. Gubser and P. J. E. Peebles, Phys. Rev. D70, 
123510 (2004). 

[16] M. Kesden and M. Kamionkowski, Phys. Rev. D74, 
083007 (2006). 

[17] M. Kesden and M. Kamionkowski, Phys. Rev. Lett. 97, 
131303 (2006). 

[18] L. Ackerman, M. R. Buckley, S. M. Carroll, and 
M. Kamionkowski, Phys. Rev. D 79, 023519 (2009). 

[19] J. L. Feng, M. Kaplinghat, H. Tu, and H.-B. Yu, J. 
Cosmol. Astropart. Phys. 07 (2009) 004 . 

[20] G. DAmico, M. Kamionkowski, and K. Sigurdson, 
arXiv:0907.1912. 

[21] J. L. Feng, arXiv: 1002.3828. 

[22] D. Tucker-Smith and N. Weiner, Phys. Rev. D64, 
043502 (2001). 

[23] F. J. Sanchez-Salcedo, Astrophys. J. 591, L107 (2003). 
[24] M. Abdelqader and F. Melia, Mon. Not. R. Astron. Soc. 

388, 1869 (2008). 
[25] A. H. G. Peter, Phys. Rev. D81, 083511 (2010). 
[26] R. A. Flores, G. R. Blumenthal, A. Dekel, and J. R. 

Primack, Nature (London) 323, 781 (1986). 
[27] R. Con, Astrophys. J. 546, L77 (2001). 



[28] V. Springel, N. Yoshida, and S. D. M. White, New As- 
tronomy 6, 79 (2001). 

[29] V. Springel, Mon. Not. R. Astron. Soc. 364, 1105 
(2005). 

[30] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 462, 563 (1996). 

[31] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 490, 493 (1997). 

[32] G. L. Bryan and M. L. Norman, Astrophys. J. 495, 80 
(1998). 

[33] J. Diemand and B. Moore, arXiv:0906.4340. 

[34] S. Kazantzidis, J. Magorrian, and B. Moore, Astrophys. 

J. 601, 37 (2004). 
[35] W. H. Press et al., Numerical Recipes in C: The Art 

of Scientific Computing (University Press, Cambridge, 

England, 2nd ed., 1992). 
[36] C. Power, J. F. Navarro, A. Jenkins, C. S. Frenk, 

S. D. M. White, V. Springel, J. Stadel, and T. Quinn, 

Mon. Not. R. Astron. Soc. 338, 14 (2003). 
[37] Y. Lu, H. J. Mo, N. Katz, and M. D. Weinberg, Mon. 

Not. R. Astron. Soc. 368, 1931 (2006). 
[38] J. S. Bullock, T. S. Kolatt, Y. Sigad, R. S. Somerville, 

A. V. Kravtsov, A. A. Klypin, J. R. Primack, and 

A. Dekel, Mon. Not. R. Astron. Soc. 321, 559 (2001). 
[39] A. F. Neto, L. Gao, P. Bett, S. Cole, J. F. Navarro, C. S. 

Frenk, S. D. M. White, V. Springel, and A. Jenkins, 
Mon. Not. R. Astron. Soc. 381, 1450 (2007). 
[40] A. V. Maccio, A. A. Dutton, F. C. van den Bosch, 

B. Moore, D. Potter, and J. Stadel, Mon. Not. R. As- 
tron. Soc. 378, 55 (2007). 

[41] A. R. Duffy, J. Schaye, S. T. Kay, and C. Dalla Vecchia, 
Mon. Not. R. Astron. Soc. 390, L64 (2008). 

[42] A. V. Maccio, A. A. Dutton, and F. C. van den Bosch, 
Mon. Not. R. Astron. Soc. 391, 1940 (2008). 

[43] D. H. Zhao, Y. P. Jing, H. J. Mo, and G. Borner, As- 
trophys. J. 707, 354 (2009). 

[44] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, M. R. 
Nolta, C. L. Bennett, B. Gold, M. Halpern, R. S. Hill, 
N. Jarosik, et al., arXiv: 1001. 4635. 

[45] D. N. Spergel, L. Verde, H. V. Peiris, E. Komatsu, 
M. R. Nolta, C. L. Bennett, M. Halpern, G. Hinshaw, 
N. Jarosik, A. Kogut, et al., Astrophys. J. Suppl. Ser. 
148, 175 (2003). 

[46] D. N. Spergel et al., Astrophys. J. Suppl. Ser. 170, 377 
(2007). 

[47] J. M. Comerford and P. Natarajan, Mon. Not. R. As- 
tron. Soc. 379, 190 (2007). 

[48] D. A. Buote, F. Gastaldello, P. J. Humphrey, L. Zappa- 
costa, J. S. Bullock, F. Brighenti, and W. G. Mathews, 
Astrophys. J. 664, 123 (2007). 

[49] R. Mandelbaum, U. Seljak, and C. M. Hirata, J. Cos- 
mol. Astropart. Phys. 08 (2008) 006. 

[50] J. M. Comerford, M. Meneghetti, M. Bartelmann, and 
M. Schirmer, Astrophys. J. 642, 39 (2006). 

[51] J. F. Hennawi, N. Dalai, P. Bode, and J. P. Ostriker, 



16 



Astrophys. J. 654, 714 (2007). 

[52] D. H. Rudd, A. R. Zentner, and A. V. Kravtsov, Astro- 
phys. J. 672, 19 (2008). 

[53] J. Annis, S. Bridle, F. J. Castander, A. E. Evrard, 
P. Fosalba, J. A. Fricman, E. Gaztanaga, B. Jain, A. V. 
Kravtsov, O. Lahav, et al., arXiv:astro-ph/0510195. 

[54] N. Kaiser, H. Aussel, B. E. Burke, H. Boesgaard, 
K. Chambers, M. R. Chun, J. N. Heasley, K. Hodapp, 

B. Hunt, R. Jedicke, et al., in Society of Photo- Optical 
Instrumentation Engineers (SPIE) Conference Series, 
edited by J. A. Tyson & S. Wolff (2002), vol. 4836 of 
Presented at the Society of Photo-Optical Instrumenta- 
tion Engineers (SPIE) Conference, pp. 154-164. 

[55] Paul A. Abell, J. Allison, S. F. Anderson, J. R. Andrew, 
J. R. P. Angel, L. Armus, D. Arnett, S. J. Asztalos, T. S. 
Axelrod, S. Bailey, et al. (LSST Science Collaboration), 
arXiv:0912.0201. 

[56] P. Predehl, G. Hasinger, H. Bohringer, U. Briel, 
H. Brunner, E. Churazov, M. Freyberg, P. Friedrich, 
E. Kendziorra, D. Lutz, et al., in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference 
Series (2006), vol. 6266 of Presented at the Society of 
Photo-Optical Instrumentation Engineers (SPIE) Con- 
ference. 

[57] A. Albrecht, G. Bernstein, R. Cahn, W. L. Freedman, 
J. Hewitt, W. Hu, J. Huth, M. Kamionkowski, E. W. 
Kolb, L. Knox, et al., arXiv:astro-ph/0609591. 

[58] E. Romano-Diaz, I. Shlosman, Y. Hoffman, and 

C. Heller, Astrophys. J. 685, L105 (2008). 

[59] S. Pedrosa, P. B. Tissera, and C. Scannapieco, Mon. 
Not. R. Astron. Soc. 402, 776 (2010). 

[60] P. B. Tissera, S. D. M. White, S. Pedrosa, and C. Scan- 
napieco, arXiv:0911.2316. 

[61] F. Governato, C. Brook, L. Mayer, A. Brooks, G. Rhee, 
J. Wadsley, P. Jonsson, B. Willman, G. Stinson, 
T. Quinn, et al., Nature (London) 463, 203 (2010). 

[62] A. R. Duffy, J. Schaye, S. T. Kay, C. Dalla Vecchia, 
R. A. Battye, and C. M. Booth, arXiv:1001.3447. 

[63] W. H. Press and P. Schechter, Astrophys. J. 187, 425 
(1974). 

[64] J. R. Bond, S. Cole, G. Efstathiou, and N. Kaiser, As- 
trophys. J. 379, 440 (1991). 

[65] S. D. M. White, G. Efstathiou, and C. S. Frenk, Mon. 
Not. R. Astron. Soc. 262, 1023 (1993). 

[66] R. K. Sheth and G. Tormen, Mon. Not. R. Astron. Soc. 
329, 61 (2002). 

[67] B. E. Robertson, A. V. Kravtsov, J. Tinker, and A. R. 
Zentner, Astrophys. J. 696, 636 (2009). 

[68] M. Oguri, K. Takahashi, H. Ohno, and K. Kotake, As- 
trophys. J. 597, 645 (2003). 

[69] A. Mantz, S. W. Allen, H. Ebeling, and D. Rapetti, 
Mon. Not. R. Astron. Soc. 387, 1179 (2008). 

[70] J. Dunkley, E. Komatsu, M. R. Nolta, D. N. Spergel, 

D. Larson, G. Hinshaw, L. Page, C. L. Bennett, B. Gold, 
N. Jarosik, et al., Astrophys. J. Suppl. Ser. 180, 306 
(2009). 

[71] J. P. Henry, A. E. Evrard, H. Hoekstra, A. Babul, and 

A. Mahdavi, Astrophys. J. 691, 1307 (2009). 
[72] A. Mantz, S. W. Allen, D. Rapetti, and H. Ebeling, 

arXiv:0909.3098. 
[73] J. Tinker, A. V. Kravtsov, A. Klypin, K. Abazajian, 

M. Warren, G. Yepes, S. Gottlober, and D. E. Holz, 

Astrophys. J. 688, 709 (2008). 
[74] J. E. Carlstrom, G. P. Holder, and E. D. Reese, Annu. 



Rev. Astron. Astrophys. 40, 643 (2002). 

[75] A. D. Hincks, V. Acquaviva, P. Ade, P. Aguirrc, 
M. Amiri, J. W. Appel, L. F. Barrientos, E. S. Bat- 
tistelli, J. R. Bond, B. Brown, et al., arXiv:0907.0461. 

[76] Z. Staniszewski, P. A. R. Ade, K. A. Aird, B. A. Benson, 
L. E. Bleem, J. E. Carlstrom, C. L. Chang, H. Clio, 
T. M. Crawford, A. T. Crites, et al., Astrophys. J. 701, 
32 (2009). 

[77] K. Ichiki, M. Oguri, and K. Takahashi, Physical Review 
Letters 93, 071302 (2004). 

[78] I. Zehavi et al., Astrophys. J. 630, 1 (2005). 

[79] L. E. Strigari, J. S. Bullock, M. Kaplinghat, J. D. Si- 
mon, M. Geha, B. Willman, and M. G. Walker, Nature 
(London) 454, 1096 (2008). 

[80] S. Vegetti, L. V. E. Koopmans, A. Bolton, T. Treu, and 
R. Gavazzi, arXiv:0910.0760. 

[81] L. V. E. Koopmans, M. Barnabe, A. Bolton, M. Bradac, 
L. Ciotti, A. Congdon, O. Czoske, S. Dye, A. Dutton, 
A. Elliasdottir, et al., arXiv:0902.3186. 

[82] P. J. Marshall, M. Auger, J. G. Bartlett, M. Bradac, 

A. Cooray, N. Dalai, G. Dobler, C. D. Fassnacht, 

B. Jain, C. R. Keeton, et al., in AGB Stars and Re- 
lated Phenomenastro2010: The Astronomy and Astro- 
physics Decadal Survey (2009), vol. 2010 of Astronomy, 
pp. 194 — h. 

[83] L. A. Moustakas, K. Abazajian, A. Benson, A. S. 

Bolton, J. S. Bullock, J. Chen, E. Cheng, D. Coe, A. B. 

Congdon, N. Dalai, et al., arXiv:0902.3219. 
[84] U. Seljak, Mon. Not. R. Astron. Soc. 318, 203 (2000). 
[85] Z. Zheng, J. L. Tinker, D. H. Weinberg, and A. A. 

Berlind, Astrophys. J. 575, 617 (2002). 
[86] K. Abazajian, Z. Zheng, I. Zehavi, D. H. Weinberg, J. A. 

Frieman, A. A. Berlind, M. R. Blanton, N. A. Bahcall, 

J. Brinkmann, D. P. Schneider, et al., Astrophys. J. 

625, 613 (2005). 
[87] Z. Zheng, A. A. Berlind, D. H. Weinberg, A. J. Benson, 

C. M. Baugh, S. Cole, R. Dave, C. S. Frenk, N. Katz, 
and C. G. Lacey, Astrophys. J. 633, 791 (2005). 

[88] Z. Zheng and D. H. Weinberg, Astrophys. J. 659, 1 

(2007) . 

[89] G. Gilmore, M. I. Wilkinson, R. F. G. Wyse, J. T. 
Kleyna, A. Koch, N. W. Evans, and E. K. Grebel, As- 
trophys. J. 663, 948 (2007). 

[90] L. E. Strigari, J. S. Bullock, and M. Kaplinghat, Astro- 
phys. J. 657, LI (2007). 

[91] L. E. Strigari, S. M. Koushiappas, J. S. Bullock, 
M. Kaplinghat, J. D. Simon, M. Geha, and B. Willman, 
Astrophys. J. 678, 614 (2008). 

[92] M. G. Walker, M. Mateo, E. W. Olszewski, 
J. Penarrubia, N. Wyn Evans, and G. Gilmore, Astro- 
phys. J. 704, 1274 (2009). 

[93] E. J. Tollerud, J. S. Bullock, L. E. Strigari, and B. Will- 
man, Astrophys. J. 688, 277 (2008). 

[94] M. I. Wilkinson and N. W. Evans, Mon. Not. R. Astron. 
Soc. 310, 645 (1999). 

[95] J. D. Simon, A. D. Bolatto, A. Leroy, L. Blitz, and E. L. 
Gates, Astrophys. J. 621, 757 (2005). 

[96] W. J. G. de Blok, F. Walter, E. Brinks, C. Trachternach, 
S.-H. Oh, and R. C. Kennicutt, Astron. J. 136, 2648 

(2008) , 0810.2100. 

[97] R. Kuzio de Naray, S. S. McGaugh, and W. J. G. de 

Blok, Astrophys. J. 676, 920 (2008). 
[98] L. V. E. Koopmans, T. Treu, A. S. Bolton, S. Buries, 

and L. A. Moustakas, Astrophys. J. 649,599 (2006). 



17 



[99] A. B. Newman, T. Treu, R. S. Ellis, D. J. Sand, [102] P. J. Humphrey, D. A. Buote, F. Gastaldello, L. Zappa- 

J. Richard, P. J. Marshall, P. Capak, and S. Miyazaki, costa, J. S. Bullock, F. Brighenti, and W. G. Mathews, 

Astrophys. J. 706, 1078 (2009). Astrophys. J. 646, 899 (2006). 

[100] J. Richard, G. P. Smith, J. Kneib, R. Ellis, A. J. R. [103] F. Gastaldello, D. A. Buote, P. J. Humphrey, L. Zappa- 
Sanderson, L. Pei, T. Targett, D. Sand, M. Swinbank, costa, J. S. Bullock, F. Brighenti, and W. G. Mathews, 
H. Dannerbauer, et al., arXiv:091 1.3302. Astrophys. J. 669, 158 (2007). 

[101] R. Mandelbaum, U. Seljak, R. J. Cool, M. Blanton, [104] R. Mandelbaum, G. van de Ven, and C. R. Keeton, Mon. 

C. M. Hirata, and J. Brinkmann, Mon. Not. R. Astron. Not. R. Astron. Soc. 398, 635 (2009). 
Soc. 372, 758 (2006). 



